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Control of charge migration in molecules by ultrashort laser pulses 
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Due to electronic many-body effects, the ionization of a molecule can trigger ultrafast electron 
dynamics appearing as a migration of the created hole charge throughout the system. Here we 
propose a scheme for control of the charge migration dynamics with a single ultrashort laser pulse. 
We demonstrate by fully ah initio calculations on a molecule containing a chromophore and an 
amine moieties that simple pulses can be used for stopping the charge-migration oscillations and 
localizing the charge on the desired site of the system. We argue that this control may be used to 
predetermine the follow-up nuclear rearrangement and thus the molecular reactivity. 

PACS numbers: 82.50.Nd, 31.15.A-, 82.53.-k 


Since the early days of quantum mechanics, the aim of 
the scientific efforts has been not only to understand the 
microscopic world, but also to use the quantum prop¬ 
erties for controlling different processes. For example, 
one may use the quantum interference and the prop¬ 
erties of the laser-matter interaction in order to exert 
control over the chemical reactivity of a molecule^ re¬ 
search field known nowadays as “femtochemistry” [l| . In 
femtochemistry, by specifically tailored femtosecond laser 
pulses, one tries to steer the motion of the nuclei in the 
molecule and thus to bring the system into the desired 
reaction channel. 

In recent years, the rapid development of the attosec- 
ond pulse techniques opened the door for studying and 
eventually controlling the electronic motion. As long 
as the notion of electronic states is valid, the potential- 
energy landscape in which the nuclei move is formed by 
the faster moving electrons. It is, therefore, very appeal¬ 
ing to use attosecond pulses that modulate the electronic 
motion such that it triggers some desired nuclear rear¬ 
rangement. Using the electron dynamics and the quan¬ 
tum coherence to induce a particular chemical process 
is the new paradigm of the emerging field of “attochem- 
istry” (see, e.g., Ref. Bi). 

This novel concept was to a large extent motivated by 
the interpretation of the results of a series of pioneer¬ 
ing experiments performed by Schlag, Weinkauf, and co¬ 
workers (see, e.g.. Refs. 0-0) showing a charge-directed 
reactivity Q in electronically excited ionic states of var¬ 
ious peptide chains. It was observed that after localized 
ionization of the chromophore site of the peptide, the pos¬ 
itive charge is transferred to the remote end of the chain 
causing a bond-breaking. Time-resolved measurement 
on a smaller, prototype molecule (2-phenylethyl-A,N- 
dimethylamine, abbreviated as PENNA) showed that the 
process takes place on the time scale of few tens of fem¬ 
toseconds B B] ■ Extensive ab initio many-body calcula¬ 
tions on PENNA suggested 0 as explanation that the 
charge-directed reactivity is concerted electron-nuclear 
dynamics: immediately after the ionization pure electron 
dynamics are triggered and the positive charge starts to 


oscillate between the chromophore and the remote amine 
end of the molecule on a few-femtosecond time scale, 
while at later times the coupling to the slower nuclear 
dynamics causes the trapping of the charge on the amine 
site and the bond breaking. 


How is the charge transferred in this pure electronic 
step? In their pioneering work till, Cederbaum and Zo- 
beley demonstrated that such a transfer, termed charge 
migration, can be solely driven by the many-body effects 
(electron correlation and electron relaxation). Due to the 
electron correlation, the removal of an electron from a 
molecular orbital will create an electronic wave packet (a 
simultaneous population of a multitude of cationic states) 
which will evolve in time. The charge migration has 
been intensively studied theoretically 12 and turned 
out to be a rich phenomenon, with many facets that are 
rather characteristic of the molecule studied (for a recent 
review, see Ref. [l7|). We note also that very recently, 
pure electronic, few-femtosecond charge oscillations were 
experimentally observed after a broadband ionization of 
phenylalanine [ll| (see also Ref. [l9|). 


In this Letter we present a scheme for controlling 
the many-body electron dynamics of the charge migra¬ 
tion process by a single ultrashort laser pulse. We 
show that by appropriately tailored infrared pulses one 
can stop the pure electronic, few-femtosecond oscilla¬ 
tion of the charge, localizing it on the desired site of 
the molecule 3-methylen-4-penten-A,A-dimethylamine 
(MePeNNA), which is a structural analogue of the men¬ 
tioned above PENNA. 


As noted already, the reason for the charge migration 
is that due to the electronic correlation, the removal of 
an electron from a particular molecular orbital populates 
several ionic states, creating in that way an electronic 
wave packet ll|. Depending on the type of the popu¬ 
lated cationic states, different mechanisms of charge mi¬ 
gration have been identified 00,0 , with the most 
common in the outer valence being the so-called hole¬ 
mixing mechanism. In the hole mixing, two (or more) 
ionic states appear to be linear combinations of two (or 
more) one-hole (Ih) configurations, representing removal 
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of an electron from a particular molecular orbital. 



FIG. 1. (color online) Ground and first excited cationic 
states of the molecule MePeNNA computed using the ab initio 
many-body Green’s function ADC(3) method [^. The next 
ionic state is located at 10.5 eV. The contributions of the 
Ih configuration (HOMO)“^ are given in red, while that of 
(HOMO—1)“^ is shown in green. The two molecular orbitals 
are also depicted. The figure is adapted from Refs. f^ . 


can be achieved experimentally without approaching the 
sudden-approximation limit of removing an electron from 
a single orbital. As we see from Fig.[Tl for preparing such 
a wave packet we need a coherent population of the first 
two ionic states only, which can be done via a laser pulse 
with photon energy centered between the two states and 
a band width sufficient to embrace both of them. Since 
the states are about 0.5 eV apart, one needs a pulse with 
duration of about 1 fs. The next ionic state is located at 
10.5 eV and therefore its population by such a pulse will 
be negligible. The initial localization site of the charge 
(chromophore or N-terminal) is determined by the rel¬ 
ative phase between the two ionization channels. The 
latter can be controlled through the ionization pulse pa- 
rameters, e.g., by chirping the pulse [2^, or by using a 
TT-pulse [24 1. 

We now pose the question whether after the electron 
dynamics are triggered and the charge starts to bounce 
back and forth between the two ends of the system, we 
can control its motion by applying a short laser pulse. 
The interaction of the system with an external electric 
field E(t) can be described (in the dipole approximation) 
by the Hamiltonian 


A show-case example of a strong hole mixing appears 
in the molecule MePeNNA. Our many-body calculations 
showed that due to the electron correlation the ground 
and the first excited ionic states of the molecule are 
a strong mixture of two Ih configurations: an elec¬ 
tron missing in the highest occupied molecular orbital 
(HOMO) and an electron missing in the HOMO—1, see 
Fig. [Hand Refs. (2^, |2l[. Therefore, if we suddenly re¬ 
move an electron either from HOMO or from HOMO—1, 
we will create an electronic wave packet, which in the 
Schrodinger representation reads (in atomic units): 

\Mt)) = ci(t)e-*^^‘|/) + C2(t)e-*^^‘| J)- (1) 

In the above expression ci{t) and C 2 (t) are the time- 
dependent (in general complex) amplitudes of the 
cationic eigenstates |/) and |J), satisfying |ci(t)p + 
|c 2 (t)P = 1 at all times. Due to the hole-mixing structure 
of the ionic states, the evolution of the system described 
by the the wave packet dm will represent an oscillation of 
the hole charge between the two involved Ih configura¬ 
tions, or between HOMO and HOMO—1. Since HOMO 
is localized on the chromophore and HOMO—1 on the 
amine group of the molecule (see Fig. |T|), the charge mi¬ 
gration will represent an oscillation of the charge between 
the two ends of the molecule with frequency determined 
by the energy difference of the two states uq = Ej — Ej, 
which is 0.55 eV, meaning that the time needed for the 
charge to reach the remote end of the molecule is only 

3.8 fs [iSlili. 

Interestingly, the ionization spectrum of MePeNNA 
suggests that an initial state of the form of Eq. CD 


H{t)=H,on-D-E{t), 


( 2 ) 


where Hion is the full electronic Hamiltonian of the ion¬ 
ized system and D = {Dx, Dy,Dz) is the vector operator 
of the dipole moment. The evolution of the system is gov¬ 
erned by the time-dependent Schrodinger equation with 
the formal solution 


|$,(t)) = Texp |$,(0)). (3) 

In practice, one usually breaks the interval of interest 
[0,t] into a large number of small increments, or time 
steps At, and Eq. is integrated numerically. 

To describe the correlated motion of 69 electrons (the 
number of electrons in MePeNNA cation) in the presence 
of a laser field is an extremely difficult problem. To solve 
it we use ab initio methods only. The cationic Hamilto¬ 
nian Hion is constructed using the so-called non-Dyson 
algebraic diagrammatic construction (ADC) scheme 221 
for representing the one-particle Green’s function. At 
the ADC(3) level, used in the present calculations, the 
Hamiltonian is consistent with the exact Green’s function 
up to the 3rd order of perturbation theory with respect 
to the Hartree-Fock (HF) reference Hamiltonian. Stan¬ 
dard DZ basis sets [23 were used for constructing the 
non-correlated reference states. 

One can include the interaction of the field by diagonal¬ 
izing the cationic Hamiltonian matrix and using the field- 
free eigenstates as a basis to expand the dipole operator, 
computing the transition dipole matrix elements. Alter¬ 
natively, when the diagonalization is very expensive, one 
can represent the dipole operator D in the many-electron 
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basis of the molecular Hamiltonian (in the present case 
this is the so-called intermediate-state-representation ba¬ 
sis [2^) and then directly propagate with the full Hamil¬ 
tonian (^ . For this purpose, or for performing the wave- 
packet propagation, Eg. ®, we used the short-iterative 
Lanczos technique 27| . Details of this technique allowing 
to study correlated-electron dynamics in systems contain¬ 
ing few 10s of electrons are given in Ref. 28|. 

Our aim is to find laser fields E{t) which can steer the 
evolution of the wave function of a system. As men¬ 
tioned above, under certain ionization conditions, the 
initial state of the MePeNNA molecule will be a linear 
combination of the two lowest eigenstates of the ionic 
Hamiltonian and, therefore, can be regarded as a two- 
level system. There are numerous protocols suggested in 
the literature for controlling quantum dynamics in two- 
level systems (see, e.g.. Refs. [iiil). Most of these 
protocols, however, are designed to perform a popula¬ 
tion inversion, that is, to invert the population distribu¬ 
tion between the two states. This is insufficient for our 
purposes, as we would like to have the freedom to choose 
the desired final distribution of the populations starting 
from any initial one. 

Very recently, we proposed a general approach for ob¬ 
taining resonant laser pulses that can drive the popula- 
tions of a two-level system in a predefined way [33l| . If we 
want that the evolution of the system follows a particu¬ 
lar quantum path, that is, that one of the populations in 
Eq. ([T]) evolves according to some control function f{t), 
i.e., |ci(t)p = f{t), then the field which can drive this 
transition takes the form 13311 : 


m = 


1 


m 


M \//W(i - /W) 


sin {ujot + If), 


( 4 ) 


where ^ is the projection of the dipole moment on the 
electric field polarization axis, and ip is the relative phase 
between the initial amplitudes, ci(0) and 02(0), of the 
populated ionic states. A convenient choice for the con¬ 
trol function is 13311 


fit) = 


1 - 
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1 -I- e“ 
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1 


1 


( 5 ) 


where Qi and a / are the initial and the final (target) pop¬ 
ulations, respectively, of one of the states and the param¬ 
eter a connects the transition time with the intensity of 
the held. 

It should be noted that the two-level model is used 
here only for obtaining the held parameters. For com¬ 
puting the evolution of the hole charge in MePeNNA 
molecule in the presence of the control pulse, the propa¬ 
gation was performed with the full Hamiltonian, Eq. ©, 
as described above. 

A convenient quantity for visualization, or for tracing 
in time and space the charge dynamics is the so-called 
hole density ll|, ll2| . The hole density is dehned as the 


difference between the electronic density of the neutral 
and that of the cation 

Q{r,t) = (4'o|p(f=')|4'o) - ($*(t)|p(r)|$i(t)), (6) 

where p{r) is the one-body electronic density operator, 
Id^o) is the ground state of the neutral, and |$i(<)) is the 
propagated cationic wave packet. The quantity Q(f,t) 
describes the density of the hole at position r and time t 
and by construction is normalized at all times t. 

Let us now examine two situations of particular inter¬ 
est for achieving control over the charge migration dy¬ 
namics, namely stopping the charge oscillations and lo¬ 
calizing the charge on one of the two molecular sites. To 
be specific, we will assume that the initial ionization of 
our test molecule, MePeNNA, is performed such that the 
electron is removed from the HOMO, i.e. the initial hole 
charge is localized on the chromophore. As discussed 
above, this will trigger pure electron dynamics in which 
the charge will oscillate between the chromophore and 
the amine site and we would like to apply a control pulse 
which will stop this oscillation and localize the charge on 
the amine group or on the chromophore. 

Such a control can be achieved by a laser pulse ob¬ 
tained via Eqs. dH) and JS]) by choosing the desired initial 
and final populations of the two states in the wave packet. 
In the case of MePeNNA, if we want to drive the system 
to a stationary state in which the charge is entirely local¬ 
ized in the HOMO, we need to choose a/ = 1, while if we 
want to localize the charge on the amine site, we need to 
take a/ = 0. The initial population is Ui = 0.4, reflecting 
the fact that an initial state with a hole localized in the 
HOMO has the form |d)i(0)) = v^OAl/) -|- J) (see 
also Fig. [T]). 

As mentioned above, the parameter a in Eq. con¬ 
trols the interplay between the duration and the inten¬ 
sity of the pulse needed to perform the transition - a 
slow transition can be achieved with a weak pulse, while 
shorter pulse will naturally need higher intensity. Since 
we want to modulate the charge migration before the 
nuclear motion will start to influence the dynamics, we 
would like to use as short pulses as possible. On the other 
hand, the high intensity of the control pulse may lead to 
undesired multiphoton processes, which can, for exam¬ 
ple, further ionize the system. Therefore, we need to bal¬ 
ance between these two factors. In the case of MePeNNA 
the minimum pulse duration for performing the needed 
transition is about 10 fs, giving an electric field strength 
which never exceeds 10® V/m. This corresponds to a 
pulse with a peak intensity of about 10^^ W/cm®, which 
is rather weak. 

The results of our full propagation accounting for the 
influence of the control field with the above parameters 
are shown in Fig. [2] The figure depicts the hole density, 
defined in Eq. ®, of the molecule MePeNNA after cre¬ 
ating the initial hole on the chromophore at time t = 0 
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FIG. 2. (color online) Time evolution of the hole density, 
Eq. di), along the molecular axis of the molecule MePeNNA, 
after an initial localized ionization of the chromophore, con¬ 
trolled with a laser pulse (shown on the right) centered at 
15 fs. The molecular axis is chosen to pass through the longest 
spatial extension of the molecule. Upper panel: the laser pulse 
is designed to achieve localization of the charge on the chro¬ 
mophore. Lower panel: the laser pulse is designed to achieve 
localization of the charge on the amine site. 


and applying a control laser pulse (also shown in the fig¬ 
ure) with the maximum of the field centered at t = 15 fs. 
Immediately after the ionization, the charge migration 
dynamics take place with the hole jumping from one end 
of the system to the other. We clearly see that between 
f ^ 10 fs and t ^ 20 fs, the time during which the system 
is exposed to the control pulse, the charge oscillations are 
nearly completely stopped and the hole becomes localized 
on the desired site - on the chromophore (upper panel 
of Fig. [2]) or on the amine group (lower panel of Fig. [IJ. 
We would like also to emphasizes that the charge stays 
put at the desired site of the molecule after the pulse is 
over, that is, the pulse is tailored such that it brings the 
system to a superposition of electronic states in which 
the density is essentially stationary. The small oscilla¬ 
tions taking place after the pulse is over reflect the fact 
that the molecule is, of course, not a perfect two-level 
system. Because of the electron correlation, in addition 
to the two mixed Ih configurations, the two involved 
ionic eigenstates contain also small contributions from 
the two-hole-one-particle (2hlp) configurations 3^ . The 
latter represent excitations on top of the ionization and 


their weight forms the missing to 1 part in the states de¬ 
picted in Fig. [TJ Through the 2hlp configurations, other 
molecular orbitals also contribute to the dynamics and 
get population, while the pulse is optimized to account 
only for the HOMO and the HOMO—1. Nevertheless, 
the suggested simple scheme for obtaining the needed 
control-pulse parameters works remarkably well in such 
a complicated system as the MePeNNA molecule. 

Let us now comment on the possible further evolution 
of the studied system. Due to its similarity with the men¬ 
tioned above PENNA molecule, one may expect a similar 
charge-directed reactivity after localized ionization of the 
chromophore. The charge starts to oscillate between the 
two ends of the molecule until the nuclear motion eventu¬ 
ally traps it at the amine site and the molecule dissociates 
by braking the bridging carbon-carbon bond. However, 
as we noted above, the nuclear motion is strongly influ¬ 
enced by the electron dynamics and, therefore, by con¬ 
trolling the charge migration we may be able to predeter¬ 
mine the nuclear rearrangement. Localizing the charge 
on the chromophore may substantially slow down or even 
prevent the dissociation, while its localization on the 
amine site will most probably speed-up the break-up of 
the molecule. 


To check this plausible hypothesis, one needs a full 
quantum treatment of the coupled electron-nuclear dy¬ 
namics. This is unfortunately currently out of reach for 
so large systems. We would like to note, however, that 
a semiclassical method that might give some hints about 
how the nuclear motion will be affected by the electron 
dynamics, was recently proposed and could be a 

good starting point for such studies. 

Before concluding, we would like to emphasize that 
the scheme presented in this work is general and not re¬ 
stricted to only stopping the charge migration oscilla¬ 
tions. Through the control function f{t) in Eq. (|T]), one 
is able to obtain the pulse shape needed to drive the sys¬ 
tem to any combination of the two electronic states, and 
thus bring the ion to the optimum initial condition for 
the desired follow-up nuclear motion. 


We hope that our study will stimulate further theoret¬ 
ical and experimental work on the possibilities to control 
chemical reactions via the manipulation of the electron 
dynamics. 
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